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Abstract 

We build Gaussian wave packets for the linear Schrodinger equation and its finite difference space semi-discretization 
and illustrate the lack of uniform dispersive properties of the numerical solutions as established in [6] . It is by now 
well known that bigrid algorithms provide filtering mechanisms allowing to recover the uniformity of the dispersive 
properties as the mesh size goes to zero. We analyze and illustrate numerically how these high frequency wave 
packets split and propagate under these bigrid filtering mechanisms, depending on how the fine grid/coarse grid 
filtering is implemented. 

Resume 

On construit des paquets d'ondes gaussiennes pour l'equation de Schrodinger lineaire continue unidimensionnelle 
ainsi que pour sa semi-discretisation en espace par differences finies. On illustre numeriquement le manque d'uni- 
formite par rapport au pas du maillage des proprietes de dispersion des solutions numeriques demontre dans 
[6]. Par ailleurs, il est bien connu que les algorithmes bi-grilles sont des mecanismes de filtrage efficaces pour 
recuperer l'uniformite des proprietes dispersives. On analyse la facon dont les solutions bi-grilles correspondant a 
plusieurs projections de la grille fine sur la grossiere, se divisent en plusieurs paquets d'ondes, chacun se propageant 
differement. On represente numeriquement ces phenomenes et on montre que ce comportement est en accord avec 
les resultats theoriques connus sur la dispersion des solutions bi-grilles. 
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Dans [6], il a ete demontre que les solutions de l'equation de Schrodingcr lineaire unidimensionnelle semi- 
discretisee en espace par differences finies sur une grille uniforme ne verifient pas les estimations dispersives 
de l'equation de Schrodingcr continue (estimations Li — L v x et gain d'une dcmi-dcrivee localement en 
espace) uniformement par rapport au pas du maillage. Ce comportement pathologique est du a l'existence 
de points critiques qui annulent l'une des deux premieres derivees de la relation de dispersion. 

Une maniere de recuperer l'uniformitc dc ccs estimations consiste a fairc un filtrage bi-grillc des donnees 
initiales. Ce mecanisme a ete introduit dans le travail pionnier de R. Glowinski [4] sur lc filtrage bi- 
grille pour le controle des ondes. Son efficacite a ete demontree rigureusement dans le contexte des 
proprietes dispersives de l'equation de Schrodingcr dans [6]. Plus precisement, tout en gardant le schema 
d'approximation semi-discret conservatif en differences finis, il s'agit de considcrer uniquement des donnees 
initiales dans un maillage quatre fois plus grossier et de les prolonger par interpolation lineaire sur la grille 
fine. 

Notre objectif principal dans cette Note est d'illustrer numeriquement les effets pathologiques des 
solutions numeriques et d'etudier dans quelle mesure ils disparaissent lorsqu'on applique les mecanismcs 
de filtrage bi-grille. 

Nous considerons des donnees initiales oscillatoires obtenues a partir d'un profil gaussien concentre 
en Fourier autour d'un certain nombre d'onde. Afin dc comparer les differentes solutions numeriques 
obtenues par le filtrage bi-grille, avant d'utiliser l'interpolation lineaire pour passer de la maillc grossicrc 
a la fine, on doit projeter les donnees initiales du maillage fin sur le grossier. On le fait de deux manieres : 
d'une part, par restriction des fonctions discretes aux points du maillage grossier et, d'autre part, en 
moyennant les valeurs correspondant aux points de la grille fine situes dans le voisinage de chaque point 
dans la grille grossiere. La premiere peut generer, entre autrcs, des solutions qui n'oscillcnt pas et nc sc 
propagent pas, concentrees en Fourier autour du nombre d'onde £ = 0, done, regulieres. La seconde peut 
donner lieu a des solutions tendant vers zero quand le pas du maillage tend vers zero. 

Dans les deux cas, on constate que ce comportement est compatible avec les proprietes dispersives 
prouvees dans [6] pour la methode de differences finies apres le filtrage bimaille. 

Problem formulation. Let us consider the 1 — d linear continuous Schrodinger equation (CSE): 

t) + d 2 x u(x, t)=0, x e K, t € M \ {0}, u(x, 0) = cp(x), x E R. (1) 

The solution of (1) is given by u(x,t) — S(t)ip(x), where S(t) = cxp(iic^) is the Schrodinger semigroup 
defined as S(t)ip(x) = (G(-,t) * <p)(x) and G(x,t) is the fundamental solution of (1). 

The solution of (f ) verifies two dispersive properties: the gain on the integrability and the local smoothing 
effect 

l|w||i«(R,i,S(R)) < c(p)|M|l2 (r) , respectively sup (— / / |<9y 2 u(x, i)| 2 dxdt) < c\\(p\\ L 2 (M) , (2) 

R K ttJ-Rj-R ' 

for (p, q) such that the admissibility conditions 2 < p < oo and 2/q= 1/2 — 1/p hold. For detailed proofs 
and higher dimensional versions of these dispersive estimates, see [8]. They play a key role in the proof 
of the well-posedness of the non-linear Schrodinger equation in L 2 (U.) or ([2], [6]). 

On an uniform grid of size h > of the real line, Qh — {xj — jh : j e Z}, we introduce the discrete 
Laplacian of the sequence ~f h = (/j)jez to be d\fj = h~ 2 (fj +1 — 2fj + fj-i), and analyze the finite 
difference semi-discrete Schrodinger equation (DSE): 

id t Uj(t) + d 2 h u 3 (t) = 0, j eZ,teR\ {0}, ^-(0) = ipj, j e Z. (3) 

Set Tlh '■= [—w/h,w/h]. The solution of (3) can be expressed by means of the inverse semi-discrete 
Fourier transform (SDFT) (cf. [6]) as Uj(t) = [cxp(itd 2 )~i^ h }j = ^ J n <p h (£) exp(itph{£)) exp(i£xj) rf£. 



2 



The symbol p h : U h -> R is defined as p h (0 = 4/i 2 sin 2 (£/i/2) and (p h is the SDFT at scale h of the 
initial data = (ipj)j e z- 

In the continuous case, the Fourier symbol of the Laplacian is = |£| 2 . Its first order derivative, 
the so-called group velocity, d^p^) = 2£, vanishes only at £ = and its second-order derivative, the 
so-called group acceleration, = 2, docs not vanish at any wave number. For the semi-discrete case, 

the following two pathologies of the symbol Ph(£) were observed (cf. [6]): 

pi. The group velocity, d^ph{£,) = 2h~ 1 sin(£/i), vanishes at £ = 0, but also at £ = ±ir/h. 

p2. The group acceleration, d^ph(£.) = 2cos(£/i), vanishes at £ = ±Tr/2h. 

The pathologies (pi) and (p2) lead to the lack of uniform dispersive properties as h — > and, more 
precisely, to the lack of local smoothing effects and of discrete Strichartz estimates, respectively (cf. [6]). 

Our main goal is to illustrate the effects of these pathologies on the numerical solutions and to analyze 
to which extent they disappear when applying the two-grid filtering mechanisms proposed in [6] and 
inspired by the pioneering work by R. Glowinski [4] on the bi-grid filtering for the control of waves. 

As proved in [6] , an efficient mechanism to recover the uniformity as h — > of the dispersive properties 
is the bi-grid algorithm introduced in [4]. This consists in solving the DSE on the fine grid of size h with 
slow initial data obtained by linear interpolation from data given on a coarser grid of size nh. We will take 
n = 2 k , with k > 1. For ratios l/2 fe , k > 2, between the two grids (the fine one of size h and the coarse 
one of size 2 k h), both pathologies (pi) and (p2) are canceled and the uniform dispersivity is recovered 
(cf. [6]). 

Define the extension operator Tfe : £ 2 (2 k hZ) — > £ 2 (h%) from the grid of size 2 k h to the one of size h, by 
linear interpolation, as follows: (T k f) 2 k j+r = (2 k -r)/2 k f 2kj +r/2 k f 2 k, J+2 k, for all j € Z and < r < 2 k -l. 

If ~f e £ 2 (2 k hZ), the Fourier transform of this extension can be written as T k f (£) = b k (£h)f 2kfl (£ i ), for 
all £ € Tlh, with weights bk(rj) — IljLi cos 2 (2 J ~ 2 ?7) vanishing quadratically at r\ = ±jn/2 k ~ 1 , for all 

1 < j < 2 k - Y . The SDFT f 1 " h is defined for £ e H 2 k h , but for simplicity, we also denote in this manner 
its extension to II ft by 7r/(2 fe_1 /i)-periodicity. 

According to [6], when k = 2, i. e. when the ratio between the two meshes is 1/4, and the numerical 
scheme is restricted to this class of filtered initial data, the dispersive properties turn out to be uniform 
as h — > 0. More precisely, there exist two constants C(p),C > independent of h such that the solution 
lt h {t) = exp(itdl)T 2 ^ h verifies 

\tf h (t)\\ LW (hz)) < C{p)\\^ h y m and sup^ fh £ \dl /2 Uj (t)\ 2 dt < C||r 2 ^|| 2 2(M) . 

_R>0 " JR lx ^ R 

Here, denotes the discrete fractional derivative, i.e. d^fj = (2n)^ 1 j Uh P^ 2 (£)/ ?l (£) exp(i^xj) d£, with 
s > 0, and the discrete £ p (ZiZ)-spaces are defined as usual (see, e.g. [6]). 

But, in practice, the initial data for the DSE are given on the fine grid of size h as an approximation of 
the initial datum of the continuous Schrodinger equation on the nodal points Xj — jh. Thus, they need, 
first, to be projected into the coarse one. We analyze two different projection operators from Qh to Q 2 k h : 




r=0 



The superscripts r and a stand for restriction and average, the two key mechanisms on which these 
projections are based. More precisely, the projection A r k restricts the function (/j)jez on the fine grid to 
those j-s that are integer multiples of 2 k . The projection A£ takes as value at the point x 2 kj an average 
of the values at the 2 k+1 — 1 points surrounding x 2 kj in the fine grid. 

These projection operators can be represented in the Fourier space, H 2 k h , as follows: 
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j=— 2*- 1 j=— 2* -1 

Behavior of solutions under quadratic dispersion relations and Gaussian initial data. Our 

aim is to describe the behavior of the numerical solutions that the bi-grid algorithm produces when the 
initial data on the fine grid is a highly concentrated Gaussian wave packet. For this, we introduce the 
Gaussian profile <7 7 (£) = y/2n/^ exp(— |£| 2 /(27)), where 

7 = j(h) such that 7/1 2 / 3 << f and 7 >> f . (6) 
Consider 770,771 < t] 2 € [— 7T,7r] independent of h and the following function depending on 770, rji, 772, 7: 

= - vo/h)x [m/ K V2 / h m (7) 

For ?7o G (— 7r, 7r], ^ represents a Gaussian profile concentrated around r/o/h truncated to [r]i/h, f]2/h] C 
fl/j. Note that even for the CSE, for an easier comparison with the DSE, we consider the same initial 
data supported in U h in the Fourier space. 

Taking into account that the dispersion relation of the CSE is quadratic, one can obtain an explicit 
representation formula for the corresponding solution uJJJ'IJ 2 . The same can be done, for more general 
dispersion relations, as it is for instance the case for the DSE, by taking the second order Taylor expansion. 
In the purely quadratic case we have: 

Qvo,h(0 = J^<lvo(Z h )> with Qvoiv) ~ Qr, (Vo) + <l'r, (Vo)(v ~ Vo) + *<4, (%)(>? - Vof ''■ (8) 

This, together with the expression of the initial data (7), allows to conclude that: 

i. Propagation: u^f} 2 propagates along the curve x(t) — x* — tq' m (r]o)/h. 

ii. Time evolution of the support: it expands like (7^ + £ 2 7(<7^' (770)) 2 ) 1 / 2 . 

iii. Time evolution of the amplitude: it behaves like c(l + t 2 7 2 (^' (r/o)) 2 )" 1 ^ 4 , where c = 1 if r/o 7^ f]i > V2 
and c = 1/2 if 770 = Vi 01 Vo = For q'^ (r)o) 7^ 0, this large time behavior agrees with the 
i~ 1/,2 -decay rate of the solutions of the CSE with Z^-data (see [6]). 

We now consider the CSE with two particular choices of initial data ip given in the Fourier space as 
follows: 

fcr (0 = T- l'n (?) + V% (0 Or (p V0 (0 = (0 . % € (-7T, 7T) . (9) 

In (9), we roughly consider initial data as in (7), but with a particular choice of 770, r?i and 772 and with 
two superposed Gaussian profiles when they are supported in any of the end points ±7r//i. The initial 
data (p n , having two picks, seems to concentrate around two wave numbers, ±ir/h, and this occurs in the 
continuous setting. But in the discrete one, since the two picks centered at ±n/h are located exactly on 
the boundary of 11^, by 27r/ft,-periodicity, they are the two halves of the same pick. The corresponding 
solutions have amplitude approximately equal to one in the physical space at t = 0. 

Denote by u„(x, t) and u Vo (x, t), 770 € (— n, n), the corresponding solutions of the CSE with initial data 
(9). Since the dispersion relation q(rj) — rj 2 for the CSE (1) is quadratic, the results above apply and the 
solution u v (x,t) splits into two blocks, uZ^(x,t) and u^(x,t), propagating along x(t) — x* ± 2irt/h, 
whereas u rt0 (x, t) propagates along x(t) — x* — 2tr) /h. 

Discrete solutions without filtering. We consider the DSE with initial data as in (9). Since the 
dispersion relation p{rj) — 4 sin 2 (77/2) is not quadratic, we cannot apply directly the above results on 
quadratic dispersions relations, but this can actually be done up to a small reminder term using the 
Taylor expansion, as mentioned before. More precisely, for all 770 G IT, we may split the dispersion 
relation as ^(77) = q^oiv) + r r)o(v)i where q Vo is the second-order Taylor polynomial about 770 and r m is 
the corresponding reminder. Consider the DSE (3) with initial data (7), denote its solution by ty~'~ 2 and 
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b) 



c) 




Figure 1. a) Initial data ip vo with r)o = w, 7r/2, 2tt/3 (blue, red, green) and their projections A£ with b) k = 1 and c) 
k = 2. In black, the corresponding weights b^. In c), the blue and red curves coincide, so that A^ip^ = A^^/a, modulo an 
exponentially small error. This also means that ip^(x4j) and W^C^-y) almost coincide. 
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Figure 2. a) Initial data ip vo with rjo = w, 7r/2, 27t/3 (blue, red, green) and their projections AJJ with b) k = 1 and 
c) fc = 2. In black, the corresponding weights b^. The amplitude of A^ip vo is smaller than the one of A. r k tp vo both in 
the physical and in the Fourier spaces, since A^ip,^ ~ fcfc(^o)A£i/9>7 and fefe(?jo) < 1- Since foi(7r) = 0, the blue curve 
in b) is almost zero. The red and green ones have amplitude b± {tt/2)\/2tz = \Z2tt/2 and &i(27r/3)\/27r = v / 27r/4. Since 
b2(n) = fe2(i"/2) = 0, the blue and red curves in c) coincide and they are almost zero. The corresponding green curve has 
amplitude fe 2 (27r/3)v^7r = V^/W. 



split it as w?nl'Jl 2 — u !^ 7 2 + u J;d,7 2 i where uJ^IJ 2 is the solution corresponding to the quadratic dispersion 
relation q Vo . It is easy to check that ||^;^ 2 ||f,2 (r)/II w rjJ,7 2 Hi 2 (R) = 0(i/Vy 3 / 2 ), which is small in finite time 
intervals iff (6) holds, so that w^if ~ u^f. 

Denote by w. K (x,t) and w m (x,t) the solutions of the DSE with initial data (9). Since p'(ir) = 0, the 
solution uv does not propagate, illustrating the lack of uniform local smoothing effect in the discrete case 
(see Figure 3 - c), the blue curve). For 770 = 7r/2, since p" (it /2) — 0, the spatial support of UW2 does not 
expand as time evolves, which agrees with the lack of uniform L'^-integrability properties (see Figure 3 
- a) and b), the blue curve). 

Bi-grid solutions. For the DSE, we now consider initial data obtained by firstly projecting the data (9) 
from the fine grid to the coarse one by one of the two projections (4) and then extending those projections 
by linear interpolation from the grid of size 2 k h to the one of size h. We identify several different cases: 

A) r)o = 7T, a = r. In the Fourier space, the projection A^ip^ has picks at 2ln/2 k h, for < \l\ < 2 k ~ 1 . 
The picks at 2lir/2 k h, 1 < \l\ < 2 k ~ 1 , are canceled as h — > since the weight vanishes exactly at those 
points. The wave packet for I — 0, being located at £ = 0, does not propagate, but it does not oscillate 
either and the discrete smoothing property holds uniformly as h — > (see Figure 3 - c), the black curve, 
for (r]o,k) = (ir,2)). This packet also decays with the rate of the CSE since both dispersion relations 
(continuous and discrete one) are tangent at £ = 0. 

5)770 € ((2r-l)7r/2 fe ,(2/* + l)7r/2 fc ), for some -2 fc - 1 + l < I* < 2 k ~ 1 -l > a = r. Set r/* := r) Q -2l*n/2 k € 
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Figure 3. Solutions of both CSE and DSE in the physical space corresponding to the initial data ip^ using the projections 
A^, a = r,a, with a) (r]o,a) = (ir/2,r), b) (rjo,ce) = (7r/2,a), c) (r)o,a) = (n, r) and d) (r]o,a) = (2-7r/3,r). Legend: green 
- solution of CSE at t = 0, magenta - solution of CSE at t = 1, blue - solution of DSE without filtering at t = 1, red - 
solution of DSE with bi-grid of ratio 1/2 (k = 1) at t = 1 and black - solution of DSE with bi-grid of ratio 1/4 (k = 2) 
at t = 1. Description: a) the red wave packets have amplitude bi(n/2) = 1/2; b) the red wave packets have amplitude 
!>?(7r/2) = 1/4 and the black one has amplitude b2(n/2) = 0; c) the red and black wave packets coincide; d) in green, the 
initial datum is asymmetric with respect to the horizontal axis. This is an effect of the projection of the continuous initial 
datum into the grid. The real part of the continuous initial datum is essentially cos(27rx/(3/i)) exp(— 7X 2 /2), with 7 = /i -1 / 4 , 
which is symmetric, indeed. But, when sampled on x = (i/i)jgz, then, for j = 0, its amplitude is 1, but, for j = 1, it is 
cos(27r/3) = — sin(7r/6) = —1/2. This explains the asymmetry of the discrete plot of the datum; the red wave packets have 
amplitude &i(2tt/3) = 1/4 and fei(-7r/3) = 3/4. The black ones have amplitude b 2 (27r/3) = 1/16, b 2 (n/6) = (6 + 3^/16, 
b 2 (-57r/6) = (6 - 3-v/3)/16 and 6 2 (-7r/3) = 3/16. 

U 2 k h . Similarly, if ij Q £ (— 7T, ~(2 k - l)n/2 k ) U ((2 k - l)7r/2 fe , ir), set <t]q := r] ± it. There are two cases: i) 
i]q = 0. Then the projection A^ip^ has picks in the Fourier space at £ = 2ln/ (2 k h), with < |/| < 2 k ~ 1 and 
the analysis follows the one in part A) (see Figure 3 - a), for (t]q ~ tt/2, k = 2)); ii) t]'q 6 (—ir/2 k , 7r/2 fe )\{0}. 
Then the projection A^(p^ has picks in the Fourier space at (r/Q + 2lir /2 k )/h, where — 2 fc_1 + 1 — s < I < 
2 k ~ x — s and s = (1 + sign(?7g))/2. The solution w r m k is a superposition of 2 k wave packets propagating 
along the lines x(t) — x* — tvi with velocity vi — 2sin(?7Q + 2lir /2 k ) /h. Observe that vi = —v_ 2 k-i + i, 
1 — s < I < 2 k ~ 1 — s, i.e. there are 2 k ~ 1 pairs of wave packets constituted by one wave packet going in each 
direction with the same velocity, and t>i_ s < v 2 k-i_ s < v 2 - s < v 2 k - 1 -i-s < • • • < v 2 k--i_ s < v 2 k-2 +1 _ s . 
If sign(?7o) 7^ sig^l^o)' then the wave packet of largest amplitude changes the direction with respect to 
the solution without filtering (see Figure 3 - d), the red curve, for (r]o,k) = (2ir/3, 1)). In this case, the 
dispersive properties arc verified uniformly. 

C) rjo = (21* + l)/2 k , for some -2 k ~ x < I* < 2 k ~ 1 - 1, a = r. Then the projection A r k ip Vo has picks 
in the Fourier space at £ = (21 + l)n/(2 k h), for all -2 fc " 1 < / < 2 fc ~ 1 - 1. The solution is a 

superposition of 2 k wave packets, each one propagating along the line x(t) = x* — tvi with velocity 
vi = 2sin((2Z — l)n/2 k )/h. With respect to B) - ii., we have vi = v 2 k-i +1 _ l for all 1 < I < 2 k ~ 2 , i.e. the 
2 fc_1 wave packets going to the left can be grouped into 2 k ~ 2 pairs constituted by two wave packets having 
the same velocity which collapse at any time t > 0, and v\ < v 2 < ■ ■ ■ < v 2 k-i. For (770, k) — (tt/2, 1), 
the two blocks in the solution do not spread as time evolves (see Figure 3 - a), the red curve), which 
confirms the fact that the bi-grid algorithm with ratio 1/2 is not enough to reestablish the uniformity of 



the dispersive estimates, as predicted by the theory in [6]. 

D) a = a, r/o € (— tt,it]. Remark that A^ip^ ~ 6fc(r;o)A^^ . Then: i) the solutions corresponding to A£ 
are of smaller amplitude than those corresponding to A£, because they involve the factor bk(f]o) < 1 (sec 
the red curves in Figure 3 - a) and b)); ii) when r] is such that bk(r)o) = 0, the projection is too 
strong, in the sense that the corresponding numerical solutions tends to zero as h — > (see Figure 2 - 
b), the blue curve for (r]o,k) — (tt, 1), or c), the blue/red curves, for (770, fc) = (n, 2), (ir/2, 2). Also the 
black curve in Figure 3 - a), for (770, k) — (ir/2, 2)). The cancellations can be explained by the fact that 
for plane waves of the form ipj = exp(ii] a j), j E Z, the identity A^T^"' 1 = holds. 

Conclusions. In this article, some subtle phenomena related to the dispersivity (and the lack of) have 
been described in the context of the Schrodingcr equation and its numerical approximation schemes. We 
have shown that, choosing appropriate high frequency wave packets, one may confirm the predictions of the 
theory (cf . [6] ) and gain new insight about the complex behavior of such solutions which, depending on the 
high-frequency around which they concentrate, may or not propagate, or propagate but not disperse. We 
also analyze and describe how these solutions arc affected by the application of bi-grid filtering techniques. 
The bi-grid solutions can be decomposed into several wave packets moving at different velocities. Some 
of them are attenuated by the bi-grid weights or by the way one projects the initial data from the fine 
grid to the coarse one. In this way, the bi-grid mechanisms may yield non-oscillatory, smooth solutions 
or solutions vanishing as h — > 0, even if the original initial data do not have those behaviors. Our 
numerical simulations confirm the predictions of the theory in [6], in the sense that the bi-gird technique, 
implemented with mesh-ratio 1/4, ensures the dispersive properties of numerical solutions. 

Open problems. There are several issues, related to those treated in this paper, worth to be explored: i) 
So far we have plotted solutions based on their Fourier description. It would be interesting to do the same, 
but using numerical approximations schemes. But then one has to take into account the changes that the 
time discretizations (cf. [6]) or some kind of boundary conditions, for example the transparent ones (cf. 
[1]), could introduce in the dispersive properties; ii) the relevance of these high frequency phenomena is 
still to be explored in the context of the non-linear Schrodingcr equations. The theoretical results in [7] 
show that dispersive schemes ensure a polynomial convergence order which improves the logarithmic one 
one gets for standard finite-difference schemes by energy methods. Also the effect of splitting methods 
([5], [3]) on these high frequency wave packets is worth investigating; iii) The issues addressed in this 
Note are totally open for non-uniform grids or Schrodinger equations in heterogeneous media. 
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